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ABSTRACT 

We describe the calibration status and data products pertaining to the GR2 and GR3 data releases 
of the Galaxy Evolution Explorer (GALEX). These releases have identical pipeline calibrations that 
are significantly improved over the GR1 data release. GALEX continues to survey the sky in the Far 
Ultraviolet (FUV, ~ 154 nm) and Near Ultraviolet (NUV, ~ 232 nm) bands, providing simultaneous 
imaging with a pair of photon counting, microchannel plate, delay line readout detectors. These 
1.25° field-of-view detectors are well-suited to ultraviolet observations because of their excellent red 
rejection and negligible background. A dithered mode of observing and photon list output pose 
complex requirements on the data processing pipeline, entangling detector calibrations and aspect 
reconstruction algorithms. Recent improvements have achieved photometric repeatability of 0.05 and 
0.03 uiab in the FUV and NUV, respectively. We have detected a long term drift of order 1% FUV and 
6% NUV over the mission. Astromctric precision is of order 0.5" RMS in both bands. In this paper 
we provide the GALEX user with a broad overview of the calibration issues likely to be confronted in 
the current release. Improvements are likely as the GALEX mission continues into an extended phase 
with a healthy instrument, no consumables, and increased opportunities for guest investigations. 
Subject headings: space vehicles: instruments — surveys — telescopes — ultraviolet: general 



1. INTRODUCTION 

We describe the calibration status and data products 
pertaining to the combined GR2 and GR3 data releases 
of the Galaxy Evolution Explorer (GALEX). GALEX is 
a NASA Small Explorer mission performing an all-sky ul- 
traviolet (UV) survey in the Far UV (FUV~ 154 nm) and 
Near UV (NUV~ 232 nm) bands with 4 - 6" resolution 
and ~ 50 cm 2 effective area. At the time of this writing, 
over 13500 square degrees of this survey are available to 
the public via the Internet from the Multimission Archive 
at Space Telescope Science Institute (MAST). 11 

The GR2 and GR3 releases have identical pipeline cali- 
brations that are significantly improved over the previous 
GR1 release. Photometric precision at a given epoch is 
0.05 and 0.03 m AB in the FUV and NUV respectively. 
We have detected a long term drift of order 1% FUV and 
6% NUV over the mission. The astromctric precision has 
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been improved to 0.5" RMS. 

The path from raw telemetry to corrected catalogs, 
images and spectra is a complicated one for the GALEX 
pipeline. This fact is illustrated by Fig. [TJ which shows 
a GALEX image in detector and sky coordinates, which 
are representative of the pre- and post-pipeline images. 
This overview highlights the data quality, errors and pe- 
culiarities that are likely to confront astronomers. As 
is true for any large dataset, knowledge of these basic 
characteristics is essential to avoid confusion. 

The GALEX mission continues into an extended phase 
with a healthy instrument, no consumables, and in- 
creased opportunities for guest investigations. While the 
survey is on-going, it is likely that we will make improve- 
ments to the calibration and discover additional sources 
of error. We expect to provide updates to this paper on 
the web at www.galex.caltech.edu, as well as in the 
literature. 

This paper is organized as follows. In Section [2] we 
summarize the design of the mission and instrument. 
In Section [3] we describe how GALEX data are recon- 
structed and what the key pipeline products are. In 
Section [4] we describe the photometric calibration of the 
data, including a discussion of saturation effects and 
drift. In Section [5] we describe the instrument resolu- 
tion across the field of view and as a function of detector 
gain. In Section [6] we describe the astrometric perfor- 
mance of the instrument. In Section we describe the 
spectroscopy calibration and present an example extrac- 
tion. Finally, in Section [8] we summarize the key points 
of the paper. 

2. MISSION AND INSTRUMENT OVERVIEW 

GALEX was launched on an Orbital Sciences Cor- 
poration (Orbital) Pegasus rocket on 2003 April 28 at 
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Fig. 1. — Left panel: A GALEX NUV "dose" image of M101 in spacecraft coordinates clearly showing the 1' dither pattern, which is 
used for all observations to smooth out the effects of the detector flat field and reduce MCP fatigue. Flight panel: The fully processed, 
background-subtracted NUV image of M101 in sky coordinates, which provides the basis for the source catalog data. 

12:00 UT from the Kennedy Space Center into a circu- 
lar, 700 km, 29° inclination orbit. This 0.5 m diameter 
instrument has a very wide 1.25° field of view sensitive to 
m AB ~ 21 in the All Sky Imaging Survey (AIS), and to 
ulab ~ 25 in the Deep Imaging Survey (DIS). The instru- 
ment provides simultaneous FUV and NUV imaging with 
a pair of photon counting, microchanncl plate, delay line 
readout detectors. These detectors are well-suited to ul- 
traviolet observations because of their excellent red rejec- 
tion and negligible background. An objective grism mode 
is also available for low resolution spectroscopy. The last 
comparable survey in this wavelength range, which ex- 
tended to niAB ~ 9 was performed by the Belgian/UK 
Ultraviolet Sky Survey Telescope on the TD1 satellite 
launch ed in 1972 and described by iBoksenberg et al.1 
(fl97l . 

GALEX makes science observations on the night side of 
each orbit (to keep sky background to a minimum) during 
"eclipses" that are typically in the range of 1500 - 1800 s. 
In the AIS, as many as 12 targets are observed during 
an eclipse, while the Medium Imaging Survey (MIS) has 
one target per eclipse and the DIS returns to the same 
target for 30 or more eclipses. As of this writing, we have 
observed over 19,000 square degrees of sky and accumu- 
lated over 4 terabytes of science telemetry. An overview 
of key performance parameters is presented in Tab. [TJ 

A GALEX detector head is shown schematically in 
Fig. [2] Individual photons incident at the cathode set off 
an isolated avalanche of current inside parallel 10-12 mi- 
cron diameter microchanncl plate (MCP) pores as they 
accelerate through a few-keV field. The resulting charge 
cloud at the rear of the MCP stack, approximately 10 7 
electrons in size, is deposited on a double-layer anode, di- 
vided, and finally measured at each of four outputs. The 
layers of the anode form a pair of orthogonal delay lines. 
Timing differences measured for the charge pulses at each 
corner of the anode are proportional to the position of 
the initiating photon event. For faint UV observations, 
the combination of low background (zero read noise) and 
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Fig. 2. — GALEX scaled tube detector head electro-mechanical 
block diagram. The NUV and FUV tube heads are nearly identi- 
cal with the principal differences being the choice of photocathode 
material (Csl for FUV and Cs2Te for NUV) and window material 
(MgF 2 for FUV, Si0 2 for NUV). Also, the cathode material is de- 
posited on the detector window in the NUV channel but directly 
on the MCP in the FUV channel. Instead of a cathode, the FUV 
window has a charged grid of wires that enhances the sensitivity of 
the detector. The NUV cathode is proximity focused on the front 
MCP, and due to its reentrant design is more prone to reflections 
from off-field sources. HVi n represents the negative applied volt- 
age, 5-6 keV, while a resistor chain (R ga p ~ 10 MOhm in series 
with the MCPs at about 30 MOhm each) sets the voltage across 
the gap (Vga P ~ 500 V) and MCP stack {Vmcp ~ 4500 V) driving 
a detector bias current Imcp ~ 50 ^A. The detector pulse outputs 
in one dimension are labeled START and STOP (or sometimes "A" 
and "C"). 



high red rejection that these detectors possess represents 
a favorable trade-off against the superior QE and field 
flatness of more conventional CCD detectors. Further- 
more, MCP detectors do not require cooling, an im- 
portant consideration in the contamination-sensitive UV 
band. In order to improve photometric performance, ob- 
servations are dithered in a spiral pattern approximately 
V in diameter (at a rate of approximately 0.5 revolutions 
per minute) that smooth out small scale variations in the 
flat field and reduce bright source fatigue on the MCPs. 

A more detail e d inst rument overview is provided in 
iMorrissev et al. | (l2005f). wh i le the mission design is 
described by iMartin et al. I (|2005t) . Further details 
of our missio n oper ation s experience are described in 
iForster et aD (|2006t ) and iMorrissevl (|2006t ). 
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TABLE 1 

Summary of measured performance parameters for GALEX. 



Item 



FUV Band 



NUV Band 



Bandwidth: a 

Effective wavelength (\ e ff): h 

Mean effective area: 

Peak effective area: 

Astrometry (la, R < 0.6°): 

Field of view: 

Photometry (ltr): 

Zero point (mOuy): 

Image resolution (FWHM): 

Spectral resolution (A/ A A): 

Spectral dispersion 

Detector background (typical): 

Total: 

Diffuse: 

Hotspots: 
Sky background (typical) : c 
Limiting magnitude (5<r): 

AIS (100 s): 

MIS (1500 s): 

DIS (30000 s): 
Linearity: 

Global (10% roll off): 

Global (50% roll off): 

Local (10% roll off):° 
Pipeline image format: 



1344 - 1786 A 
1538.6 A 
19.6 cm 2 
36.8 cm 2 at 1480 
±0.59" 
1.27° 
±0.05mxs 
18.82 
4.2" 
200 

1.64A-arcsec _1 

78 c-s" 1 
0.66 c-s~ 1 -cm -2 
47 c-s" 1 
1000 c-s- 1 

19.9 
22.6 
24.8 



1771 - 2831 A 
2315.7 A 
33.6 cm 2 
61.7 cm 2 at 2200 
±0.49" 
1.25° 

±0.0377^3 

20.08 
5.3" 
118 

4.04A-arcsec _1 

193 c-s" 1 
1.82 c-s-'-cm -2 
107 c-s- 1 
10000 c-s 1 

20.8 
22.7 
24.4 



18000 c-s" 1 
91000 c-s- 1 
114 c-s" 1 303 c-s" 1 

3840 X 3840 elements with 1.5" pixels 



a Includes wavelengths with effective area 

value. 1 * Wavelength-weighted: A e // — "zLL 

,-2 „,.-l X -l 



601 photons-s 
point sources. 

3. IMAGE RECONSTRUCTION OVERVIEW 



at least 10% of the peak 
These correspond to 246 and 
respectively. ^ These are worst-case values for 



The GALEX detector output is a list of time-ordered 
40-bit photon positions and pulse heights. These data are 
delivered to the instrument Data Processing Unit (DPU), 
which adds a timing stamp every 5 ms (or every 20 fj,s 
in uncompressed mode) . The time stamp is derived from 
the spacecraft clock, which has a drift of the order of 
1 s each week that is corrected with ground updates. 
Data are stored in the 2.8 GB ring buffer of a Solid State 
Recorder (SSR), which is down- linked to a ground station 
approximately four times each day. 

The ground pipeline reconstructs a list of photon posi- 
tions and arrival times into an image. Broadly speaking, 
the photon lists are corrected for detector effects and im- 
age motion, projected into sky coordinates, binned into 
images, flat fielded, and finally background subtracted. 
These final image products ar e then scanned by the open- 
source SExtractor program (|Bertin and Arnouts II1995I ) 
to produce catalogs of source positions and fluxes. The 
steps along this path are described in more detail in the 
following sections. 

3.1. Static Calibration 

The static calibrations are those used to assemble raw 
detector event positions into accurate locations in the 
GALEX focal plane. They also flag events that are not 
photon generated, such as the periodic "STIM" pulses 
that are always electronically superimposed (in parallel 
with real photon events) on the detector anodes, and 
regions on the MCPs that are dominated by hotspots. 

Each photon event in the GALEX science telemetry 
is a heavily encoded group of fields that, when decoded 



TABLE 2 

Key data contained in the GALEX RAW6 format 

PHOTON LISTS. PHOTONS FROM THE DETECTOR ITSELF ARE 
NOT TIME STAMPED BUT RATHER TIME ORDERED, AND THE 
TIME STAMP IS ADDED BY THE DPU BEFORE BEING 
WRITTEN TO THE SSR. 



Column 


Bits 


Description 


t 


40 


Event detection time a 


XAmC 


12 


Raw X detector position in FEE pixels 


YAmC 


12 


Raw Y detector position in FEE pixels 


X A 


5 


Wiggle (the phase of the TAC) 


X b 


3 


X axis coarse clock 


Yd 


3 


Y axis coarse clock 


Q 


5 


Detected pulse height 



a In high resolution mode, the timing precision uses 40 bits to 
attain 20 /^s resolution, however this mode is rarely used. For 
most data, only 32 bits arc significant, attaining 5 ms resolu- 
tion. 

and assembled, pro duces key values liste d in Tab. [2] In 
the GALEX design (jJelinskv et al Ir2003l ) , a free-running 
coarse clock subdivides the large detector areas into an 
8x8 grid. Detector event pulses are measured asyn- 
chronously while this clock runs, so that each position 
(proportional to the measured arrival time difference of 
pulses from either side of the delay line anode on each 
axis) comprises an integer number of coarse clock steps 
plus a fine position (essentially the coarse clock phase). 

In order to construct a raw source position (JT, Y) from 
these data, one would apply the following formulas: 



X — XAmC 

Y = Y AmC - 



- olXb 
(3Y B , 



(1) 
(2) 
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where the constants a and (3 are in the range of 2000 
(solved by minimizing source widths across the field of 
view), Xb and Yg are the coarse clock counters, and 
XAmC and YAmC are the "fine" position measurements 
representing the coarse clock phase. 12 After the detector 
positions are assembled, the following static corrections 
are applied (in order) to correct the raw positions to a 
linear representation of the actual source position in the 
focal plane: 

• Centering and scaling: Coarsely scales the raw 
positions from bins to microns. The scale factors 
arc different for each axis. 

• Wiggle: The GALEX coarse-bit sub-array ap- 
proach allows degenerate combinations of coarse 
and fine components for a given detector location. 
Each of these may have a slightly different elec- 
tronic non-linearity correction owing to the dif- 
fering phases of the time-to-amplitude converter 
(TAC) for each event. The image of a single point 
source on the detector is thus blurred by the super- 
position of all of the TAC data, an effect referred to 
as wiggle. This contrasts with smaller, more con- 
ventional time-to-amplitude conversion systems. 
These do not utilize a coarse counter bit and trans- 
late all electronic non-linearities into source posi- 
tion errors rather than source widths. With care- 
ful part selection, "wiggle" in the flight units has 
been minimized considerably. In order to elimi- 
nate it completely, each detected event includes 
a 5-bit value Xa, which represents where in the 
TAC range each measurement was made. Resid- 
ual electronic non-linearities are then removed us- 
ing lookup tables for combinations of (X,Xa) and 
(Y, Xa)- Both axes use the same Xa data because 
the X and Y pulses are correlated in time for a 
given detector position. 

• Walk: This transformation corrects position-vs- 
pulsc-hcight (Q) errors. A lookup table correction 
applied in a similar manner as wiggle, reflective 
of the fact that the electronics can not perfectly 
compensate for pulse sizes that vary by a factor of 
10 across the detector. 

• Spatial non-linearity (distortion): A lookup 
table corrects high voltage detector edge effects, 
which cause electrons to bend from their optimal 
path, distorting the image especially near the edge 
of the field. 

• Hotspot masking: Each detector has numerous 
small, high count rate features resulting from a 
combination of microscopic MCP debris and chan- 
nel defects that combined generate about half of 
the intrinsic detector background. These spots 
come and go, and are best handled by direct mask- 
ing, which eliminates only about a percent of the 
detector area. New masks arc generated roughly 
annually; hotspots are generally easy to pick out 
of images because they "reverse dither" (since they 

12 The notation derives from the ends of the delay line "A" and 
"C" , or sometimes "START" and "STOP" in other work, thus the 
position is proportional to a timing difference of A — C or "AmC." 



are fixed in detector space) and create ~ 1' circular 
artifacts in GALEX sky-space images. In the data 
pipeline, the hot spot masks create localized re- 
gions of zero effective sensitivity in detector space, 
one reason why dithered observations arc preferred 
for these detectors. Depending on pointing accu- 
racy in a given observation, the dither typically 
maps these regions into "donuts" of reduced (but 
not zero) response on the sky. 

3.2. Dynamic Aspect Correction 

The fully-corrected event positions in spacecraft co- 
ordinates accurately record where each photon landed 
in the GALEX focal plane. The next step is to recon- 
struct the spacecraft pointing on the sky as a function 
of time in order to remove the dither pattern from the 
observed data. This is the 3-parameter "aspect," which 
includes the right ascension, declination, and roll of the 
telescope boresight. For imaging mode data, a precise 
estimate is arrived at by refining the spacecraft Atti- 
tude Control System (ACS) pointing estimate with NUV 
data 13 binned into 1 s slices. These snapshot images are 
each matched to the ACT catalog (jUrban et al. 1 1 19981 ) 
of known stars, 14 and the average error vector from the 
matches to the catalog stars are used to establish a cor- 
rection to apply to the ACS estimate. The source po- 
sitions are then transformed to sky coordinates (J2000) 
using a Gnomonic projection based on the pointing at 
the interpolated time. 

Grism data undergoes a similar process, but the soft- 
ware can not match spectra to known stars at this stage 
of processing. A relative aspect solution is performed 
based solely on the instantaneous spectral image posi- 
tions, and astrometric refinement is added at a later 
stage. 

Once sky coordinates have been assigned to each pho- 
ton, an "extended photon list" is computed for internal 
use by the pipeline in forming the final data products. 
In cases where detailed, time-dependent behavior of one 
or more objects over the course of an exposure is of in- 
terest, the extended photon list can be written out to a 
so called "x-file," by special request. 

3.3. Pipeline Products 

Images and catalogs from various stages of pipeline 
processing reveal many characteristics of the complex, 
multidimensional GALEX dataset and are included with 
each data release as pipeline products. A description of 
key products is provided here for reference. 

The corrected photon data are binned into images and 
flux calibrated in order to form the basis for the GALEX 
catalogs. Count images (-cnt), which contain aspect- 
corrected detector counts per sky pixel are scaled by 

13 The process will run with FUV-only data, although owing to 
the generally decreased numbers and brightnesses of stars in FUV 
data, some parameters such as the interpolation interval need to 
be adjusted. 

14 Occasionally, matches to known stars cannot be found, usually 
due to pointing errors. In these rare cases, a good Point Spread 
Function (PSF) is achieved by performing the aspect refinement 
relative to any sources of emission detected in the field that arc 
reasonably compact. This is called relative aspect refinement mode 
and results in images with good PSFs but poor absolute astromc- 
try. Roll is not refined in relative mode. 
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the effective area and exposure time to form the flux- 
calibrated intensity (-mi) images. This step is accom- 
plished with the relative response (-rr) image, which con- 
tains the effective exposure for each pixel: 

rr(X, Y) = T,f =0 flat(X(t), Y(t)) x (1 - f dead (t)) * dt(t), 
where flat(X(t),Y(t)) is the detector flat field 15 
shifted to properly align with the dithered image time 
slice, fdead is the global electronic dead time based on 
the measured STIM pulser rate during each time slice, 16 
and dt(t) is the fraction of each 1 s time slice during which 
data was collected (i.e. HV was on and the detector was 
otherwise in a nominal observing mode). Building the 
effective exposure time in this way properly handles the 
fact that the edges of the image receive less exposure 
than the middle, and that the sources each sample a 1' 
region on the flat field. The high resolution relative re- 
sponse image {-rrhr) 1 which properly matches the count 
and intensity map pixel scales, is up-sampled from the 
6" relative response array. The intensity image is simply 
the count image divided by the high resolution relative 
response image. 

As a final step, the background in each intensity im- 
age must be estimated and removed. Since the GALEX 
detectors have negligible intrinsic background and the 
UV sky is dark, estimation of the true background is 
not properly handled by the code within SExtractor. In- 
stead, we create a custom threshold image for SExtractor 
that allows a bypass of the normal process for distin- 
guishing low rate sources from background. The typi- 
cal sky background count rates in high Galactic latitude 
GALEX pointings are ~ 10" 4 (~ 10" 3 ) photons-s -1 - 
arcsec -2 in the FUV (NUV), thus it is necessary to use 
the full Poisson distribution. The background is deter- 
mined in a set of large bins, 192" on a side. In images 
with exposure times greater than 10,000 s, the bins are 
96" on a side. For the Poisson distribution, the proba- 
bility of observing greater than or equal to k events for 
a mean rat e x is given by th e incomplete gamma func- 
tion Pfc (x) (|Press et al. II1992I) . The average is computed 
itcratively within each background bin by clipping out 
high pixels where Pk(x) < 1.35 x 10 -3 , a probability 
equivalent to 3cr for a Gaussian distribution. Even af- 
ter clipping, bright or extended sources can still bias the 
background determination, therefore, a 5x5 median filter 
is applied to the background bins. The resulting back- 
ground values are linearly interpolated to the full resolu- 
tion of the images. We have determined that the back- 
ground may be biased by a few percent due to contamina- 
tion from the wings of objects. To address this problem, 
background estimation and source extraction software is 
run twice on each GALEX image. The second run uses 
source position data from the first iteration to mask out 
source pixels and recompute the sky background and 
background-subtracted intensity images. The resulting 
background image {-skybg) is subtracted from the inten- 
sity image to form the background-subtracted intensity 

15 The flat field has a value of 1 in regions for which a zero-point- 
magnitude source generates exactly 1 deadtime-corrected photon- 
s _1 . 

16 Since the electronic STIM pulses occur at a precisely known 
rate, the fraction that are missing from the event list over a given 
interval will be known. This is the fractional dead time, which 
corrects the global rates only and docs not account for local MCP 
saturation. 



image (-intbgsub), which is the image scanned by SEx- 
tractor. 

There remain issues with the background at the 5% 
level that affect the MIS survey at the faint end. Fur- 
thermore, there are numerous artifacts not removed by 
the current pipeline that the user should be alert for. 
These include: bright star glints in the NUV images 
near the edge of the detector window; large out-of-focus 
pupil images resulting from reflections in the imaging 
window and dichroic beam splitter; a 1' diameter skirt 
around bright NUV sources (at the few-percent level in 
total energy) resulting from the proximity-focused pho- 
tocathodc; and photoemission and diffraction from the 
quantum-efficiency (QE) enhancing grid wires on the 
FUV detector window that appear as linear features 
around bright sources. Work is being done to flag and 
remove these features automatically in the pipeline soft- 
ware, but at present they amount to a small fraction of 
the survey coverage area (~ 1%). 

The primary pipeline image products and descriptions 
are listed in Tab. [H while catalog products are listed 
in Tab. 03 All of the catalog output products are gener- 
ated by SExtractor, with the exception of the very useful 
merged catalog (-meat), which combines FUV and NUV 
data in one table and adds derived columns including 
the detector coordinates of a source, detector intrinsic 
pulse height, number and brightness of nearby neighbors 
to a source, and many more. Matches between the bands 
are made by assigning a match probability to each pair 
of sources using the source separations and positional 
uncertainties. Sources with signal-to- noise (S/N) ratio 
of less than 2 or a separation greater than 3" are not 
matched across bands (but do appear in the merged cat- 
alog). GALEX fields are uncrowded enough that the 
probability of a mis-match is very small, however a mul- 
tiple candidate chain technique is used to assure that 
candidates are always matched the same way regardless 
of the order in which the process is started. 

The GALEX data products are made available to the 
general public through the GALEX Science Database 
at MAST. A dedicated website provides access to the 
pipeline output files via a web-based browser, as well as 
other derived products, such as JPEG images. The cat- 
alog data are also stored in a commercial SQL Server 
database engine fine-tuned for astronomy applications, 
which provides sophisticated query interfaces support- 
ing, for example, efficient spatial and color searches. The 
GALEX MAST database also contains enhancements 
over the photometric pipeline outputs, e.g. spatial cross- 
matching to the Sloan Digital Sky Survey (SDSS) cata- 
log. Also, since adjacent GALEX fields overlap slightly, 
a primary resolution algorithm deals with the problem of 
unique sample selection. The formal GALEX footprint 
definitions of the GALEX surveys are published in the 
National Virtual Obse rvatory's online Footprint Service 
(|Budavari et al. Il2006f) . On this site, 17 one can not only 
search and visualize the GALEX sky coverage but also 
compute the exact area in square degrees or intersect it 
with other surveys. 

4. PHOTOMETRY 
4.1. GALEX Magnitude Reference System 

17 http:/ /voservices. net/footprint 
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TABLE 3 

Binary tables generated by the GALEX pipeline. 



Extension Description 

-cat Primary SExtractor source catalogs for each band 

-feat SExtractor catalog of NUV extractions based on FUV coordinates 

-gauss Listing of the Gaussian filter elements used by SExtractor for image convolution 

-meat Merged catalog for both NUV and FUV extractions 

-neat SExtractor catalog of FUV extractions based on NUV coordinates 

-sexcols Listing of SExtractor catalog column names 

-scxparams Listing of SExtractor parameters used 



The GALEX photometric calibration is tied to the 
Hubble Space Telescope (HST) through reference spec- 
tra sup-rjhed_jrithe_CALSPEC database and described 
bv lBohlin et al. I ()2001l ). Almost all of the white dwarfs 
in this database are too bright for a photon counting 
instrument such as GALEX in direct imaging mode. As 
such, we are restricted to the dimmest sources in the cat- 
alog, which are all ultimately based on the International 
Ultraviolet Observer (IUE) calibration. All of these data 
have been re processed for th e HST calibration program, 
as detailed in lBohlinl (|1996P l. The GALEX calibration is 
thus traceable to a hydrogen white d warf model sca led 
using the ground based photometry of (|Landolt II1992D in 
combination with several white dwarf transfer standards 
measured with the HST Faint Object Spectrograph. The 
combination of models and measurements was used to de- 
fine corrections to the IUE calibration, thus transferring 
the HST calibration and allowing a fairly large group of 
standards from the IUE library to enter the CALSPEC 
database. Of these, we have observed 18 as detailed in 
Tab.|4j however we use one of the dimmest, LDS749b, as 
our primary standard for direct imaging. For the most 
part this white dwarf is in the linear range of the GALEX 
detectors - bright enough for high S/N measurements in 
just a few minutes, but dim enough not to saturate the 
microchannel plates. At the edge of the field in the NUV, 
where the detector gain is lowest, we apply a ~ 20% cor- 
rection for local saturation. This effect was discovered 
relatively recently as our database of white dwarf obser- 
vations has expanded to include a much finer sampling of 
the detectors in regions where the gain is an issue. 18 At 
the time of this writing we are planning observations of 
the (dimmest) CALSPEC white dwarf, LB227, in order 
to improve our confidence in the NUV zero point. Unfor- 
tunately this source is not much dimmer than LDS749b, 
but since we believe the saturation of LDS749b is mild it 
should provide very useful verification of our saturation 
correction. 



4.2. Magnitude Zero Point 

18 Since the detector MCPs are clamped at the edges, there is 
less chance for a charge packet at the edge of the field to spread 
laterally while traversing the stack. At the detector center, where 
there are small gaps between the plates, charge can spread across 
the pores. Thus the edges of the detectors tend to saturate more 
quickly than the centers. This effect can be seen for bright objects 
with multiple measurements by plotting the observed magnitude 
against the unsaturated detector pulse height, #uv_Q, tabulated 
in the merged catalog for each source. This quantity is the normal 
unsaturated pulse height of the detector at the location the source 
was measured. Note that the effect is minimized for large measur- 
ing apertures because the core of the PSF is primarily affected. 



GALEX use s the AB magnitude system of 
lOke and Gunnl (fl983h with the FUV and NUV 
magnitudes defined as follows: 

m uv = m0 uv - 2.5 log /[/v 

Here, fuy is the dcad-time-corrected count rate for a 
given source divided by the flat field map, which has a 
median value of order 0.85, and mOuv is the zero point 
that corresponds to the AB magnitude of a 1 count-s _1 
(cps) flat-ficld-corrected detection. There are indepen- 
dent relative response functions, and the zero points are: 

m0*w = 18.82 (3) 
m0 NUV = 20.08. (4) 

The GALEX project has maintained these zero points 
constant; they are based on the imaging- mode bandpass 
determined during ground calibration. Improvements in 
calibration have been incorporated by adjusting the flat 
field, which is a stacked sky background image modified 
by an interpolated grid of white dwarf measurements. 
Since GALEX has a large field of view, this has the ef- 
fect of moving the positions in the field of view for which 
sources with zero point magnitudes correspond to a true 
1 count-s _1 at the detector. For the most part this is not 
a major consequence since the flat field peak-to-valley 
variations for most of the useful field amount to about a 
factor of 2 even in FUV, and just 15% in NUV. There is 
thus only a modest effect on the estimation of the true 
expected noise for a given source with field position. In 
order to estimate any systematics resulting from applica- 
tion of the new flat field, a large sample of sources com- 
mon to both the GR1 and GR2/GR3 releases have been 
compared. We find that on average sources in GR2/GR3 
are shifted 0.04 and 0.01 m AB fainter in FUV and NUV, 
respectively, compared to GR1. One must of course be 
careful applying these figures to any individual measure- 
ment because the flat has changed shape as well as its 
average value. 

We have measured a drift in GR2 / GR3 photometry of 
0.003 ±0.001 mag-year" 1 FUV and 0.016 ±0.0005 mag- 
year -1 NUV. This was identified with two independent 
means. In the first method (yielding the quoted results), 
we compared the GR2/GR3 photometry of 197 FUV and 
547 NUV deep field sources between magnitudes 15 and 
20 measured at approximately yearly intervals over the 
course of the mission. We have checked this result against 
the LDS749b white dwarf data, however this second re- 
sult required considerable reprocessing. By design, the 
GR2/GR3 sample of white dwarf measurements through 
2005 shows no drift. Each measurement was given equal 
footing in order to improve the spatial sampling, which 
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TABLE 4 

HST White Dwarf Standards observed by GALEX. 



Star 


m FU v a 


m N uv a 





(2000) b 


<5(2000) b 


GD50 


11.98 


12.57 


3" 


48 m 


50.4 s 


-1° 


T 


26.8" 


HZ4 


14.53 


14.50 


3" 


55 m 


22.0 s 


It'- 


47' 


18.2" 


HZ2 


12.86 


13.25 


A h 
1 


12 m 


43.6 s 


ll 


51' 


49.9" 


G191B2B 




10.17 


5" 


5 m 


30.7 s 


52° 


if)' 


49.3" 


GD108 


! 9 


12.77 


10 h 


gm 


47.3 s 


-8° 


26' 


28.5" 


HZ21 


12.55 


13.13 


I2 h 


13 m 


56.3 s 


32° 


56' 


30.8" 


GD153 


11.33 


11.91 


I2 h 


57 m 


2.4 s 


22° 


1' 


53.6" 


HZ43 


10.75 


11.36 


13' 1 


16 m 


21.8 s 


29° 


5' 


57.8" 


HZ44 


10.02 


10.27 


13 h 


23 m 


35.0 s 


36° 


r 


58.8" 


GRWp70d5824 


12.14 


12.75 


ld> h 


38 m 


49.2 s 


70° 


17' 


4.6" 


SAP041C 




16.69 


14 h 


51 m 


58 s 


71'- 


13' 


17.4" 


BD+33°2642 


10.51 


10.47 


15 h 


51 m 


59.9 s 


32° 


56' 


54.8" 


SAP177D 




18.35 


15'' 


59 m 


13.5 s 


47° 


36' 


41.4" 


SAP330E 




17.84 


16' 1 


31™ 


33.8 s 


30° 


8' 


46.3" 


LDS749B C 


15.57 


14.71 


21 h 


32 m 


16.4 s 


0' J 


15' 


15.0" 


G93-48 


12.14 


12.39 


21 h 


52 m 


25.2 s 


2° 


23' 


17.9" 


NGC7293 


10.93 


11.70 


22 h 


29 m 


38.5 s 


-21° 


9' 


47.9" 


LTT9491 


16.09 


14.58 


23 h 


19 m 


35.4 s 


-18° 


51' 


32.0" 



a Magnitudes shown arc predictions based on the GALEX bandpass and 
publicly available reference spectra from the HST CALSPEC database at 
http: / / www.stsci.edu/instrumcnts / observatory / cds / calspcc.html b Coordi- 
nates are GALEX-mcasurcd in the NUV. One exception is BD+33°2642. 
which was only observed in grism mode; the coordinate provided is from 
IBohlin et al. I J2001D . C This white dwarf is the primary GALEX standard. 



only included 20 measurements per year early in the mis- 
sion. In 2005, we identified the NUV white dwarf satura- 
tion effect in low gain areas, and also recovered the FUV 
detector from a major anomaly. In light of these cir- 
cumstances, we implemented an extensive recalibration 
campaign. Over 1500 NUV measurements were added, 
but only 137 in the FUV because the detector was still 
completing its recovery. Thus the GR2/GR3 flat, gener- 
ated with all of the white dwarf samples from the mis- 
sion through 2005, is biased to 2005 in both bands and 
very strongly so in NUV. However, areas of the flat field 
sampled only in, say, 2003, will reflect the proper cali- 
bration for that epoch, which means that there is an ad- 
ditional systematic built into the flats for isolated areas. 
To address the drift, we reprocessed all of the GALEX 
LDS749b white dwarf data using a new flat field based 
on a grid of 165 LDS749b field measurements made over 
2 days in 2006. We used a nearby field star of mod- 
erate NUV count rate rather than the bright LDS749b 
to generate the NUV flat (for FUV we used the white 
dwarf), mitigating concerns about the NUV saturation. 
The results of the white dwarf and deep field analyses are 
consistent with each other and are presented together in 
Fig. El The data are consistent with a linear drift (al- 
though less certain in the FUV band). We define the 
time of zero correction to be at eclipse 12000 based on 
the weight of that epoch in generating the GR2/GR3 
flats. 

We have not included an FUV LDS749b data point 
for 2005 in Fig. [3] This is because of a second issue we 
discovered during this analysis, which is that the 2005 
FUV white dwarf data reads bright by about 5% when 
reprocessed with the 2006 single-epoch flat field. This 
is a result of the way the measurements were made in 
2005 while the FUV detector was completing its recovery 
from an anomaly. During that period (approximately a 
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Fig. 3. — The measured photometric drift for deep field sources 
between magnitudes 15 and 20 over the course of the GALEX mis- 
sion. FUV data is shown in blue: diamonds for the deep field 
sources and circles for LDS749b. NUV data is shown in red: tri- 
angles for the deep field sources and circles for LDS749b. The 
deep field data shown represent the mean binned in 1000 eclipse 
segments. The drift is defined to be zero at eclipse 12000, which 
occurred in summer 2005 during the extensive white dwarf observ- 
ing campaign upon which the current calibration is mostly based. 
It is not clear if the trend is due to contamination in the instru- 
ment, or due to a mild global gain decrease (fatigue) in the detector 
microchannel plates. 

month), the detector was being cycled on and off for 
only very short observation periods (~minutes) while its 
stability was being tested. Since the detector was ramped 
to high voltage literally minutes before the observation 
began, it was still warming up during the observations, 
offsetting the calibration slightly. This means that there 
is a field-dependent systematic built into the GR2/GR3 
FUV flat field of approximately 5% in the sense that 
most (about 77%) objects are actually brighter in the 
FUV than they appear. NUV is unaffected. 
We have not determined if the observed drift is due 
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Fig. 4. — GALEX photometric "curve of growth," which shows 
the median magnitude offsets relative to the largest aperture used 
for the white dwarf LDS749b. NUV measurements have been re- 
stricted to high pulse height areas of the detector to minimize any 
effect due to saturation. 

to contamination or a global MCP fatigue. Since the 
performance of each MCP pore is sensitive to its total 
exposure, We would expect fatigue to affect NUV more 
since it typically detects photons at about 10 times the 
FUV rate. Furthermore, contamination would normally 
be expected to preferentially affect the FUV band. On 
the other hand, the effect we observe appears to affect 
bright and dim objects the same. Whichever effect is at 
play, we expect to apply a correction automatically in 
the pipeline during a future calibration update, however 
it is not accounted for in GR2/GR3. 

4.3. Relative precision 

For accurate and repeatable measurements of unre- 
solved sources, aperture magnitudes arc usually best 
since they are fixed from measurement to measurement. 
One needs to take care to choose an appropriate aperture 
for the task at hand. Smaller apertures are more useful 
for unresolved sources, but larger ones become desirable, 
even for point sources, when comparing the photometry 
of objects measured at the center and edge of the field 
(where the PSF is degraded), or for bright objects that 
have saturated cores. For the purposes of this paper, 
the SExtractor MAG_APER_4 magnitude is used unless 
otherwise noted. This 12" diameter aperture is a good 
compromise between the very large apertures required 
to gather almost all of the GALEX flux in the extended 
PSF (as is done for the white dwarf calibration) and the 
small ones that are optimal for minimizing the contribu- 
tion of background. The relationship between the various 
standard SExtractor aperture magnitudes and the large 
tailor-made apertures used for white dwarf calibration 
are presented in Fig. 2] 

In the GR2/GR3 release, the scatter in white dwarf 
measurements has been reduced to about 0.02 mag RMS 
in NUV and 0.045 mag RMS in FUV. This compares to 
a scatter of order 0.1 mag in both bands for the GR1 re- 
lease. We attribute the improvement in precision to the 
way the flat field has been generated, using white dwarf 
measurements to modify the flat locally rather than us- 
ing the average to scale a sky background field. As with 
any instrument, care should be taken when interpreting 
the magnitudes of objects expected to have significantly 
biased spectral content compared to the white dwarf. 



These may have more scatter due to optical effects across 
the field of view. 

We have also verified the performance of the flat field 
by means independent of the partially saturated (in the 
NUV) white dwarf standard. For this purpose, we have 
examined the sample of data from the DIS described 
in Section 14.21 This comprises a series of fields with 
several dozen full-cclipsc-length observations, many at 
different roll angles (thus sampling multiple detector lo- 
cations with the same source). By matching compact 
(< 10" FWHM) sources with no neighbors closer than 
12" in the co-added catalogs to the individual visit cat- 
alogs, a visit database of sources and magnitude errors 
(with the co-added catalog considered to be the truth 
catalog) has been generated. A comparison of high and 
low gain detector measurements for a given source re- 
veals the count rate required to cause local saturation. 
For both detectors, the onset of saturation is at approxi- 
mately uiab ~ 15 (corresponding to 34 counts-s -1 FUV 
and 108 counts-s" 1 NUV) in the lowest gain regions of 
the detector. We thus restrict further analysis in this 
Section to the range 15 < tuab < 17. 

With the unsaturated magnitude range determined 
over the whole field, we can proceed to compute the 
photometric repeatability as shown in Fig. [5] Here, we 
present histograms of the 12"-diameter aperture magni- 
tude errors in each band. The Poisson error associated 
with each detection is of order 0.5% to 1.5%, while the 
total scatter we measure in the bright deep field sources 
is ±0.05 m A B RMS in FUV and ±0.03 m A B RMS in 
NUV. These figures agree well with the white dwarf re- 
sults, thus concerns about the (mild) LDS749b NUV sat- 
uration correction are reduced. If the white dwarf sat- 
uration had not been properly accounted for, the flat 
field would have been adversely affected in the NUV low 
gain areas, which would in turn increase the photomet- 
ric scatter of dimmer (definitely unsaturated) deep field 
sources in those same regions, and thus we would not 
expect good agreement in the scatter between the white 
dwarf and deep field samples. 

The deep field visit magnitude errors are also presented 
as a function of the co- added "truth" magnitude in Fig.[6l 
which serves to illustrate an interesting bias effect on the 
faint end of the 1500 s depth MIS. Here, the predicted 3a 
errors (including the flat-ficld-limitcd values derived for 
bright sources) are overplotted, showing good agreement 
with most of the observed distribution but failing on the 
faint end. This effect is similar in magnitude but oppo- 
site in direction for each detector. In the FUV, sources 
are detected in the co-added image that are not detected 
in the visit images except when the Poisson errors tilt the 
detection to the bright side (the Eddington bias), thus 
faint FUV detections are biased bright by about 0.1 m^s 
compared to the "truth" catalog. Similarly in the NUV, 
more sources are detected in the co-added image, how- 
ever the effect of the undetected sources is to force an 
overestimate of the NUV background of about 5% (with 
large scatter). In this case, faint NUV detections are bi- 
ased dimmer by about 0.1 vhab compared to the "truth" 
catalog. In each panel, solid overplotted lines show the 
median of the magnitude errors in each magnitude bin, 
while dot-dash lines show the same quantity after the 
visit magnitudes have been re-computed using the co-add 
background estimation. This correction is indistinguish- 
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clm KUV = +/- 0.050 for 412 measurements. 



-0.20 -0.10 0.00 0.10 0.30 

dm (visit-coadd) 

Fig. 5. — GALEX photometric precision as determined by re- 
peated measurements of bright but unsaturated sources in the 
GALEX deep fields. These measurements are in good agreement 
with the results determined directly from repeat measurements of 
the white dwarf calibrators. 

able from the solid line median in the FUV panel because 
the background is negligible, but completely cancels the 
faint end bias in the NUV. Since the effects of these two 
biases are in opposite directions for each band, they have 
a more significant effect on the UV color, about 0.2 irab 
at the faint end of the distribution. 

For completeness, we present a similar analysis of the 
photometric errors in the 100 s AIS sample relative to 
cross-matched sources in overlapping MIS fields. As il- 
lustrated in Fig. [TJ we find that the AIS magnitudes are 
biased bright at the faint end relative to the MIS for 
both detectors, and that corrections to the AIS magni- 
tudes using MIS backgrounds have no effect at all. This 
result stems from the shallower depth of the AIS survey, 
which is source rather than background-limited. 

4.4. Linearity 

There are two sources of photometric non-linearity in 
the instrument: global non-linearity resulting from the 
finite period required for the electronics to assemble pho- 
ton lists (during which time new events are locked out), 
and local non-linearity resulting from the MCP-limited 
current supply to small regions around a bright sources. 
The local effect complicates the NUV white dwarf cali- 
bration measurements. 

Global effects arc easily measured using the on-board 
"STIM" pulser, which electronically stimulates each de- 
tector anode with a steady, low rate stream of electronic 
pulses. These are resolved as point sources in the cor- 
ners of the dose (-scdose) images. Since the real rate of 
STIM pulses is accurately known, the measured rate is 
used by the pipeline to scale the effective exposure. The 
count rates for all sources in the field arc simultaneously 
corrected. This effect is typically about 10% in NUV 
and a few percent in FUV, however it can become quite 
significant (~ 50%) for the brightest fields. 



Local non-linearity can not be corrected well with the 
existing calibration. It affects not only the measured 
count rate of an individual bright source but also the 
source shape and position. Furthermore, the effects vary 
around the field of view with the gain of the detector 
(as exemplified by the white dwarf saturation correction 
described in Section 14. ip . We have used standard stars 
to estimate the local non-linearity in each band as shown 
in Fig. [H This data may be used to approximately cor- 
rect bright sources that are within the measured range. 
Note that the measurements presented in the Figure are 
for 34.5" (aper_7) and 3' diameter apertures; the ef- 
fect will be much worse when small apertures are used. 
The NUV detector is more robust to bright sources be- 
cause its photoelectrons are proximity-focused and thus 
present a larger source image (with lower count density) 
to the MCP. Local photometric non-linearity is strongly 
source-size dependent. The saturation shown in Fig. [8] 
is a worst-case scenario (using stars). Wide scatter in 
measurements of bright sources is an indication of the 
variation of the onset of saturation at different detector 
locations. 

We have increased the local count rate limit to 
30,000 cps per source in NUV, which corresponds to ap- 
proximately 9th magnitude. This new limit allows most 
of the sky to be observed. The FUV local limit is main- 
tained at 5,000 counts-s^ 1 , which is not currently a sur- 
vey constraint. Global detector limits are set in hardware 
(rather than the planning system) to 100,000 counts-s _1 , 
a limitation based largely on the boundaries of ground 
testing and engineering judgment about safe rates of gas 
generation inside the sealed detector tubes. Areas such 
as the galactic plane and the Magellanic Clouds are still 
restricted by the 100,000 cps global limit 

4.5. The GALEX Deep Fields 

The specific case of faint DIS photometry presents its 
own special challenges. The relatively large GALEX 
PSF affects the phot ometry of faint s ources, mostly due 
to source confusion (jXu et al. I I2005D . As an alterna- 
tive to the GALEX pipeline for faint source extrac- 
tion ! _a_F3ayesiaii_arjproach has been attempted by oth- 
ers (jGuillaume et al. Il2006[ ) to estimate the UV flux us- 
ing prior information from external, well-resolved images 
(mostly from visible wavelengths). The basic idea is to 
generate an image with prior sources convolved with the 
GALEX PSF, and to estimate what scaling factors to 
apply that will best fit the UV image. The implemented 
method is based on a maximum likelihood parametric 
optimization in the presence of Poisson noise. Scaling 
factors were extracted via an iterative procedure to solve 
the set of non-linear equations based on the Expectation- 
Maximization (EM) method. In order to test the relia- 
bility of the photometry extracted using this method, 
simulated objects (with a large range of UV flux) were 
added to the GALEX DIS images. For an NUV DIS 
image with 60 ksec exposure, it was possible to recover 
the photometry of simulated sources with an accuracy of 
0.3 m.AB to NUV~25.5, extending the faint-end limit of 
the survey by about a magnitude. 

5. RESOLUTION 

There are numerous contributors to the width of the 
GALEX PSF, including the optics, detectors, and ground 
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Fig. 6. — GALEX photometric repeatability as a function of magnitude as determined by repeated measurements of sources in the 
GALEX deep fields (with individual measurements of order 1000 s). In each panel, heavy dashed lines are the 3<r values that would be 
expected given the flat field precision and Poisson counting statistics, while the X symbols show the actual value for each magnitude bin. 
Solid overplotted lines show the median of the magnitude errors in each magnitude bin, while dot-dash lines show the same quantity after 
the visit magnitudes have been re-computed using the co-add background estimation. 
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Fig. 7. — GALEX AIS photometric repeatability as a function of magnitude as determined by repeated measurements of sources in 
overlapping MIS fields. In each panel, heavy dashed lines are the 3cr values that would be expected given the flat field precision and Poisson 
counting statistics, while the X symbols show the actual value for each magnitude bin. Solid overplotted lines show the median of the 
magnitude errors in each magnitude bin, while dot-dash lines show the same quantity after the visit magnitudes have been re-computed 
using the co-add background estimation (in the AIS this correction does not have a significant effect). At the very faint end, there is an 
apparent bias toward visit magnitudes brighter than co-added magnitudes due to Eddington bias in the source-limited sample. 



pipeline. The reconstruction algorithms for GR2/GR3 
have improved to the point that the instrumental con- 
tributions are now the largest share, particularly in the 
proximity-focused NUV channel. As we will show, the 
GALEX detectors dominate in this area and have signif- 
icantly position-dependent performance. 

To verify the fundamental instrument performance 
from on-orbit data, some bright stars were analyzed in- 
dividually outside the pipeline. These results show per- 
formance that is consistent with or better than what was 
measured during ground tests. We have also verified the 
end-to-end performance including the pipeline by stack- 



ing images of stars from the MIS survey that were ob- 
served at different locations on the detector. The results 
of these composites are shown in Figures |H] and [TUJ Per- 
formance is reasonably uniform except at the edge of the 
field, where it is significantly degraded. Some remaining 
degraded areas of the images may yet be improved with 
updates to the static calibration. Note that DIS surveys 
co-add data with many roll angles, so the DIS PSF is 
considerably more symmetrical across the field of view. 

Many of the regions with degraded resolution correlate 
with low detector gain, an effect illustrated in Fig. [TTJ 
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Fig. 8. — Measured count rate versus predicted count rate for the 
white dwarf standards measured to date by GALEX using 34.5" 
(aper_7) and 3' diameter apertures. These measurements repre- 
sent a worst case because they are for point sources. The non- 
linearity will be much worse when measuring fluxes using small 
apertures, since this effect is concentrated at the core of each im- 
age. 







FUV 


MIS P3F 


Composite 








i»l-3n,ai) 


aew-i&saj 


ten {-4, as) 


SaB(E,3S) 


ST9(14^a] 








• 


• 


# 


• 


• 








nj-e.a 




■BU!(E.BQ,T.lB] 


■»Hll!(4.M,T4J 1) 








lsf-SH.au) 






M0!-4,ao) 


S40(e,3a) 


4C9<io>ae] 






% 


• 




• 


I 


• 


# 








Trai!(4.3^Miq t 












a»(-si P nj 


B4W(-.34,1B] 


7»(- 1B>14} 


W0!-4,W) 


H0(e,i4) 


4C0<io>ie] 


5B9(3D,1H] 




• 


• 


• 




• 


# 


# 


• 


71H1I:(5. 19,4,7 4J 


71Hll:(4.71,4Ja] 

q-iaw 


7WBH:(4.n^4r4DJ 7 

q- lm 


KHU:(4. 7^444] 


■«BH:(4.»AV3] 


ntHU:(5.EB^E4] 


ntHU:(5. 39^73] 
flj-5.1 


niBH:(5.aiA30} 

»-7.0' ' 


100,-31,0) 




405!- 1M) 


ewi-4*J 




«70(1W) 


070(90,4 


495(33,4] 


• 


• 


• ' 


• 


# 


# 


• 





T1HII:(4.4<V4J 1} 

q-13.1 


L4.4 


TfBH:(4.CA4SQ T 
q-L4S 


■BH:(4.7M^!Q 
q-*B 


1BH:(4.4«ASQ 


TfHH:(C.B0,4J4] 
1~fl.B 


ITOKfo a 1,4.78] 
q-5.1 


nmil:(B.ai^4SQ] 
0J-7.B 


3m!-»,-3) 


Tin (-a4, -4) 


600!- 14,-4) 


-4,-4) 


ei^B,— ti 


«0(lfc-O) 


440(34,-4) 


M0(9J,-4) 


# 






# 


• 


# 


• 




mH:(4. 90,4^1] 




r*BH:(4.iA4.iQ] f 


iBH:(4.aa^4.iT} 


1BH:(4. 14,4a 4] 




nrBH:(4.7S^4JJT} 


mH:(B. 00^4.7 H) 


q-ns 


q- lt.i 


q-LU 


q-u» 




q-ifti 


q-ifti 


H-5.B 


Jl«(-Wv-13) 


B»(-34.-14} 


e»(-J4.-i4} 


M01-5.-14) 


740(5- L4) 


5W(J0.-lfl 


7Wt».-14} 


3W(!R-lfl 


• 


• 


• 


• 


# 




• 


• 




imm !(4, aSrWJB) 














0-14* 


q-!KS 


q-w« 


q-L0>B 


q-L4* 


q-ni 


q-u>4 


q-5.i 




4»(-ai,-aa3 


e»(-14,-33J 


O60( -4,-35) 


M0(5,-M) 


549(JC>-i3B} 


449(34,-^3 3} 


1*135,- If) 




• 


# 


♦ 












71H1I 7*4J3 1} 




IHII:(1.7T^7} 


1HII:l!4,3T,4Jje} 


rs2H;t5.tf*4.ae 


nraii:(5.i4^CB} 


nraii:(5.B4#jJ4} 




q-i7S 


q-i&s 


q-iu 


q- UV.T 


q-is* 


q-ioa 


B>-7.7 






B»(-l»>-B13 


S!0!-4,-ai) 


310(5,-31) 


lBfl(Jft-50} 














• 


• 










T1BII:(i.H^.7B} T 


KH»:f4.#»,4^S 


TOMlh.Sl^i!! 


T1BII:(5.»#^3] 










H-ia-t 




q-L4* 


q- Li* 







Fig. 9.— Composite PSFs for the GALEX FUV channel in 
detector coordinates. The individual images are exaggerated for 
clarity. Upper data values in the format N@(X,Y), where N is the 
number of stacked images and (X,Y) is the position in minutes of 
arc from the field center. The FWHM is specified in seconds of arc 
for each detector axis. Q is the typical pulse height at the specified 
detector location, in bins. The pulse height scale is not extremely 
well calibrated, but it is linear with 10 bins ~ 1 X 10 7 electrons). 

The relative detector gain (or charge, 'Q') in the region of 
measurement for each source is tabulated in the merged 
catalog. This gain measurement is obtained from a look- 
up table that represents the detector in an unsaturated 
condition. The actual gain measured during the obser- 
vation is available in the pulse height (-scq) image data 
product. 

It is also evident that the GALEX PSF has structure a 
significant distance from the center. We have character- 
ized this by building a composite PSF of the white dwarf 
calibrator LDS749b as shown in Fig. [TU The effects of 
local non-linearity are also illustrated here by the highly 



Fig. 10.— Composite PSFs for the GALEX NUV channel in 
detector coordinates. The individual images are exaggerated for 
clarity. Upper data values are N@(X,Y), where N is the number 
of stacked images and (X,Y) is the position in minutes of arc from 
the field center. The FWHM is specified in seconds of arc for 
each detector axis. Q is the typical pulse height at the specified 
detector location, in bins. The pulse height scale is not extremely 
well calibrated, but it is linear with 10 bins ~ 1 X 10 7 electrons). 
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Fig. 11. — The full width at half maximum of the GALEX 
PSF based on composite images for each detector as a function of 
intrinsic detector pulse height. In general, low pulse height regions 
are near the edge of the detector. 

saturated white dwarf HZ 44, the brightest we have ob- 
served with GALEX in imaging mode. The effect of the 
saturation is strongest in the core of the PSF, causing 
the wings to be strongly emphasized. 

6. ASTROMETRY 

The GALEX GR2/GR3 astrometric precision repre- 
sents an improvement of order 20% over GR1 (and much 
more compared to the Early Data Release). The dynamic 
pipeline transformation is generally not a limiting factor 
as it had been in earlier releases. The static transfor- 
mation (primarily the distortion map) is also performing 
very well for most of the field of view. Improvements 
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Fig. 12. — GALEX FUV and NUV composite encircled en- 
ergy curves for two white dwarf standards. A composite of 
LDS749b measurements forms the unsaturated (solid) curve, while 
the brightest standard that we have observed, HZ 44, defines the 
saturated curve. The reference count rates for the saturated curves 
are the predicted values; measured rates are much lower as indi- 
cated by the low encircled energy values. 



could still be made that would not only improve the 
global astrometric precision but also the point spread 
function (PSF) at the edge of the field. 

The distortion maps for the GR2 / GR3 calibration are 
the result of several flight-based iterations on the pre- 
ffight calibration. The first, which was imple mented for 
the G R1 data release, used Tycho-2 catalog (|Hog et al. I 
l2000f ) reference stars to refine the distortion map for each 
detector independently. An iterative technique was used 
until the measured GALEX positions converged with the 
reference positions, a technique that relied on the as- 
sumption that the distortion across a ~ 1' (~ 1 mm) re- 
gion of the detector (the footprint of the dither pattern) 
is reasonably represented by a single average value. This 
technique worked well, however there were relatively few 
FUV stars available to perform the analysis (a few thou- 
sand) . The FUV-NUV relative astrometry retained some 
noticeable (~ 1") misalignments. We have addressed this 
issue with cross-matched SDSS dat a. Since the SPSS as- 
trometry is better than 0.1" RMS (jPier et al. 1120031 ). it 
can be regarded as essentially perfect for our purposes. 
The sample we gathered contains 69642 objects spectro- 
scopically classified as quasars (SDSS SPECClass = 3, 
a convenient point source sample compiled for other ap- 
plications). Objects were considered matches if they fell 
within 2.5" of each other. There was no discrimination 
based on location in the GALEX field of view. We binned 
these results into a relatively coarse grid (~ 1'), the in- 
tent of which was to "stretch" the existing plate solution 
to remove low spatial frequency errors. The result of this 
improvement is that the GR2/GR3 absolute astrometric 
precision achieved by the GALEX MIS within 0.6° of the 
field center for bright (signal-to-noise ratio > 10) point 
sources is 0.49" RMS in the NUV and 0.59" RMS in 
the FUV, with similar relative precision between the two 
bands. 

In Fig. [13] we present histograms of the position errors 
for high S/N SDSS quasars detected by GALEX. The 
most noticeable image defect resulting from the current 
astrometric solution is the degraded PSF at the edge of 



Relative to SDSS for GALEX 
S/N>10 end field-of-view 

FUV AIS, o=0 

NUV AIS. tr=0. 



", N=26091 
N = 6202 
, N = 21029 




Fig. 13. — Distribution of absolute astrometric errors in the 
GALEX GR2/GR3 survey data for quasars with S/N> 10 in each 
band matched against the SDSS. 



Relative to SDSS for GALEX sources with S/N>10 

FUV MIS: o=0.55 M . N= 4777, R<0.5° 

FUV MIS: <7=0.60 M . N= 1425, 0.5°<R<0.6 

FUV MIS: a=0.93", N= 202, R>0.6° 




Fig. 14.— FUV GR2/GR3 MIS absolute astrometric perfor- 
mance for field center and field edge distributions of SDSS quasars 
with GALEX S/N> 10 




Fig. 15.— NUV GR2/GR3 MIS absolute astrometric perfor- 
mance for field center and field edge distributions of SDSS quasars 
with GALEX S/N> 10 

the field. As shown in Figures [14] and HH the astro- 
metric precision in the outer 5' of the detector, while 
degraded, still performs nearly at the level of the central 
30'. Generally speaking, matches between GALEX and 
other surveys can be safely performed with a 2" window 
as indicated by Fig. [TBI 
The most promising future improvement to the astro- 
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Fig. 16. — The percentage of GALEX-SDSS matches as a func- 
tion of matching radius. An MIS sample with S/N> 10 and field 
radius less than 0.6° was used for this analysis. Note that the 
matched sample is limited to 2.5" maximum error, although based 
on the shape of the distribution this does not appear to be an issue. 

metric calibration probably involves a new pipeline mod- 
ule that would scale the amplitude of the distortion maps 
by a factor determined from measurements of the STIM 
pulses at the corners of each image. Much of the neces- 
sary data to eliminate this "breathing" mode is already 
available. Another more complex possibility involves re- 
processing the reference star data in ~ 20 s intervals to 
reduce the footprint of each source sample on the detec- 
tor and thus more finely sample the detector grid. Such 
an analysis would be most useful near the edge of the 
detector where the distortion map is changing rapidly. 

7. SPECTROSCOPY 

The spectroscopic mode utilizes a CaF2 grism (a prism 
with a ruled surface) that can be moved into the converg- 
ing beam of the telescope to form simultaneous (slitless) 
spectra of all sources in the field in both bands. Many 
spectra will normally overlap, a problem addressed by 
rotating the grism on its axis to different position an- 
gles during repeat observations of the same field. Each 
spectrum has an "undeviated wavelength point" around 
which it will rotate as the grism position angle is changed 
(thus any given spectrum will have an optimal set of 
grism position angles depending on the location of nearby 
sources). As the position angle changes, the undeviated 
wavelength point remains fixed in the field as illustrated 
by Fig. [TTJ The resulting spectra can be combined to 
eliminate confusion. The standard GALEX spectroscopy 
pipeline is designed for isolated point sources, which 
are detected in the GALEX direct image observations. 
Extended or overlapping sources, or sources which do 
not appear in the direct image catalog require a special 
pipeline reduction with modified parameters. 

A relative position correction as a function of time is 
applied to the spectral image data using a technique very 
similar to that used for direct image photons. Using the 
brightest sources from the catalog, an aspect correction 
is applied to the entire exposure by fitting four param- 
eters: shift in dispersion direction, shift in spatial di- 
rection, spacecraft roll angle, and grism position angle. 
A spectrum points North (blue to red) when the grism 
position angle is 0° and towards the East when it is 90° . 

Bright sources may show grism orders ranging from -4 
to +5 across the field, however only two orders in each 




N 

Fig. 17. — A composite image of GALEX grism spectra collected 
with multiple grism position angles. As the grism is rotated in its 
mount, each spectrum spins around its own undeviated wavelength 
point to form a "pinwheel" in the composite, illustrating the op- 
timal grism angle (producing an unconfused spectrum) for each 
source. 

band are typically extracted. For NUV, the highest re- 
sponse is in first order with a total system peak response 
of 40 cm 2 . For FUV, the strongest order is 2nd with 
a peak response of 20 cm 2 . The dispersion angles jfuv 
and jnuv in arcseconds from the undeviated wavelength 
point for FUV 2nd order and NUV 1st order are given 
by: 

~<FUV = -4530.7 + 7.3S59A - 3.9254 X 10~ 3 A 2 + 7.4750 X 10~ 7 A 3 
~fNUV = -882.1 + 7.936 X 10 _1 A - 2.038 X 10~ 4 A 2 + 2.456 X 10~ 8 A 3 

where A is the wavelength in Angstroms, accurate to a 
few tenths of an arcsecond. 

Using the catalog positions, image strips (aligned to 
the dispersion direction) are created for each source us- 
ing the corrected photon positions. Fig. [IS] shows how 
a portion of the field image translates into a dispersion- 
aligned image strip for each band. The Figure shows 
image strips covering 1700 arcseconds in the dispersion 
direction to illustrate the weaker orders. Only 900 arc- 
seconds (from -100 to +800 arcseconds relative to the 
undeviated point) covering the two strongest orders in 
each band are stored during actual data reductions. The 
width of each image strip is 78 arcseconds. Fig. [18] also 
shows the undeviated point in each image strip (marked 
with an x) which occurs at about 1700 A in 1st order in 
both bands. This is the point where the spectral images 
for all position angles intersect. 

The spectral shape, strength, and position of neighbor- 
ing sources are estimated from the direct image catalog 
and arc masked from each image strip. Additional mask- 
ing is done on each strip when the exposures are coad- 
ded to eliminate other artifacts not predicted by known 
neighboring sources, which may be due to reflections, 
ghosts, hot spots, or unidentified sources near the edge of 
the field. The response map is then applied and the back- 
ground is fit and subtracted at each wavelength point. 
An aperture extraction (summation) is computed for the 





Fig. 18. — Example bright object GALEX spectra. Circular images at top show the reconstructed FUV and NUV spectral image data. A 
rectangular region of interest identifies the same spectrum in each band, while an X marks the undeviated wavelength point around which 
the spectra would spin if the grism were rotated. The rectangular regions are enlarged and annotated in the middle of the frame, where 
wavelength ranges and spectral orders are identified. In the bottom panel, the calibrated spectrum is shown in blue for the FUV portion 
and red for the NUV portion. 



primary spectrum using the brightest order in each band 
(NUV 1st and FUV 2nd). An optimal extraction is also 
computed using a PSF model. A 20 Angstrom region is 
set to zero in the spectrum at around 1820 A indicating 
the separation between FUV and NUV. 

A list of sources for which spectra will be extracted is 
determined by estimating what the final S/N of a given 
source will be after coadding all grism observations. Any 
source which will achieve a S/N per resolution element 
in the coadded result of > 2 in FUV or > 3 in NUV is 
selected for spectral extraction. For multi-exposure data 
reductions, the masked image-strip photon and response 
data are coadded for each source. For both individual 
and coadded spectra, a response correction is applied as 
a function of wavelength and the spectra are rebinncd to 
a linear scale. The primary (NUV 1st order and FUV 2nd 
order), secondary (NUV 2nd order and FUV 3rd order), 
and optimal (the primary spectrum with PSF weighting 
in the spatial direction) extractions are stored in 3 ar- 



TABLE 5 
Typical signal-to-noise ratio 
for unresolved galex 
spectra at mis depth 
(1500 s). 





FUV S/N 


NUV S/N 


16 


10.4 


21.2 


18 


3.4 


6.1 


20 


0.79 


1.22 


22 


0.14 


0.20 



rays in the final reduced data file. The estimated S/N 
for unresolved sources in single-orbit spectra is shown in 
Tab. [3 

The usable wavelength range is approximately 1300 to 
1820 A for FUV (2nd order) and 1820 to 3000 A for NUV 
(1st order). The average dispersion is 1.64 and 4.04 A 
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per arcsecond for FUV (2nd order) and NUV (1st order), 
respectively. The average resolution (A/AA) for a point 
source (assuming a 5 arcsecond PSF) is about 200 (8A) 
and 118 (20A) for FUV (2nd order) and NUV (1st order), 
respectively. 

Although there is evidence that the response may vary 
with position on the detector or with the grism position 
angle, a single wavelength-dependent calibration is used 
for each band. These variations may introduce a calibra- 
tion error of about 3% in the body of each spectral order 
and about 10% near the edges of each order. 

8. CONCLUSIONS 

We have described the GALEX performance and cali- 
bration results for the GR2 and GR3 data releases. These 
releases have identical pipeline calibrations that are sig- 
nificantly improved over the previous GR1 release. The 
analysis we have presented provides an overview of the 
data quality, errors and peculiarities that are likely to 
confront astronomers. 

The relative photometric precision is 0.05 and 0.03 
mAB in the FUV and NUV respectively. We have iden- 
tified several photometry issues that will be improved in 
future calibrations. The first is a systematic drift with 
time of order 1% in FUV and 6% in NUV, total, over the 
period of the released data (~4 years). The sense of the 
drift is that objects appear to be growing fainter with 
time. The second issue is that the current FUV calibra- 
tion relies on data collected in a non-standard mode dur- 
ing the 2005 FUV recovery, and this data appears bright 
by 5% compared to other measurements of the same ob- 
jects collected in the standard mode. The third issue 
is that the NUV zero point is defined by a white dwarf 
calibrator that is at least partially saturated. While we 



have evidence that the saturation is mild and correctable, 
we are planning confidence-building observations of an 
additional white dwarf, the dimmest in the CALSPEC 
catalog, to help resolve the issue. 

The astromctric precision has been improved, about 
20%, to 0.5" RMS. There is evidence that it may be im- 
proved further by accounting for the detector "breath- 
ing" mode already measured with existing data from an 
external pulser that runs in parallel during all detector 
operations. 

All of our performance parameters continue to meet 
flight requirements. As the mission continues, we nat- 
urally expect the calibration to evolve. The instrument 
itself continues to function nominally in spite of a se- 
ries of flight anomalies especially with the FUV detec- 
tor. We have modified operations and flight software 
over the mission to respond to these issues. While some 
of them have significantly affected the FUV observation 
efficiency, none give any indication of affecting the instru- 
ment performance. We believe there is room to increase 
the count rate limits for both detectors (especially FUV), 
and to the extent possible we plan to explore higher rate 
fields closer to the Galactic plane after the primary sur- 
veys are completed. 



GALEX (Galaxy Evolution Explorer) is a NASA Small 
Explorer, launched in April 2003. We gratefully acknowl- 
edge NASA's support for construction, operation, and 
science analysis for the GALEX mission, developed in 
cooperation with the Centre National d'Etudes Spatialcs 
of France and the Korean Ministry of Science and Tech- 
nology. 
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TABLE 6 

A SUMMARY OF GALEX PIPELINE PRODUCT ARRAY CONTENT AND FORMAT. 



Product 


Extension 


Format 


Units 


Scale 




Description 


Artifact Flags 


-flags 


480 X 


480 




12"-pixel" 


-l 


Artifact flags for each pixel 


Background-Subtracted Intensity 


-intbgsub 


3840 x 


3840 


counts-s -1 


1.5"-pixel 


-l 


Intensity-Sky Background 


Count 


-cnt 


3840 x 


3840 


counts 


1.5"-pixcl 


-l 


Aspect-corrected image, no fiat field. 


Dose 


-scdose 


2250 x 


2250 


counts 


3"-pixel" 


i 


Detector space, corrected for walk, wiggle and distortion. 


Exposure 


-exp 


960 x 


960 


s 


6"-pixel~ 


i 


Exposure, uncorrected for flat field or dead time 


Intensity 


-int 


3840 x 


3840 


counts-s -1 


1.5"-pixel 


-l 


Count image divided by relative response image 


Movie 


-movie 


480 x 


480 


counts 


12"-pixel" 


-l 


Count map slices at 16 s intervals 


Object Mask 


-objmask 


3840 x 


3840 




1.5"-pixcl 


-l 


Areas masked during sky background estimation 


Pulse height 


-scq 


2250 x 


2250 


bins 


3"-pixel" 


i 


Average pulse height per pixel 


Relative response 


-rrhr 


3840 x 


3840 


s 


1.5"-pixel 


-i 


Effective exposure, corrected for dither, flat field and dead time 


Relative response 


-rr 


960 x 


960 


s 


6"-pixel~ 


i 


Effective exposure, corrected for dither, flat field and dead time 


Sky Background 


-skybg 


3840 x 


3840 


counts-s -1 


1.5"-pixel 


-i 


Sky background estimate 


Threshold 


-wt 


3840 x 


3840 




1.5"-pixel 


-i 


Background-subtracted intensity map divided by estimated Poisson noise 



